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Abstract: We present two new aspects of Extensive Air Shower (EAS) development universality allowing to make 
accurate estimation of muon and electromagnetic (EM) shower contents in two independent ways. In the first case, to get 
muon (or EM) signal in water Cherenkov detectors it is enough to know the vertical depth of shower maximum and the 
total signal. In the second case, the EM signal can be calculated from the primary particle energy and the zenith angle. In 
both cases the parameterizations of muon and EM signals are almost independent on primary particle nature, energy and 
zenith angle. 
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Introduction 

Mass composition of ultra-high-energy cosmic rays 
(UHECR) can be studied only indirectly with large EAS ar- 
rays. The contemporary measurement of longitudinal and 
lateral shower characteristics in hybrid experiments like 
the Pierre Auger Observatory [1] provides the possibility 
to combine several primary mass sensitive EAS parame- 
ters (such as depth of shower maximum and muon shower 
content) to achieve the best primary particle mass discrim- 
ination. Unfortunately, the lack of reliable information 
on hadronic interaction properties at these energies causes 
large uncertainties in the simulations of EAS characteris- 
tics and in turn brings large uncertainties in mass composi- 
tion analysis results (see e.g. recent review [2]). 

In this paper we propose two simple, independent and accu- 
rate methods to determine muon and EM shower contents 
in hybrid experiments and briefly discuss a possible way to 
test and adjust interaction models properties in a primary 
mass independent way. We also hope that the proposed 
EAS-universality-based correction of the interaction mod- 
els will allow to perform mass composition analysis with 
the use muon EAS content in less interaction model depen- 
dent manner. 

The present study is performed making use of 
around 50000 showers, generated with COR- 
SIKA 6.735 [3]/QGSJET II [4]/FIuka [5] (see [6] for 
full list of references and more details) and COR- 
SIKA 6.900/EPOS 1.99 [7]/FIuka for R-^ spectrum in 



the energy range lg(_E/cV) = 18.5 — 20.0 and uniformly 
distributed in cos^ 9 in zenith angle interval 9 = Q° — 65°. 
EM component thinning was set to 10~^, the observation 
level was at 870 g/cm^, geomagnetic field was set to the 
value of the site of the Auger Observatory in Malargtie. 
The expected signal S in Cherenkov Auger-like detectors 
was calculated according to the sampling procedure 
described in [8, 9] with the use of the same GEANT 4 
lookup tables as in [9]. Differently from [9] in this work 
the muon signal S'p includes only signal from muons 
crossing the Cherenkov detector, while signal from EM 
particles, originating from muon decays, is included in the 
EM signal. 

1 Showers at the same XJ^^^ 

Of all aspects of universality of shower development, we 
will be interested only in dependence of EM and muon sig- 
nals on the distance of shower maximum to the ground 
and on the zenith angle. Let's consider the Auger-like 
experimental setup [1] and ground-plane signal in water 
Cherenkov detectors at 1000 meters from the shower core. 

Comparing shower characteristics dependence on the ver- 
tical depth of shower maximum ^,^iax oi^^ finds a very in- 
teresting property. Clearly, to have the same ^,^jax ^^ aver- 
age proton shower has to be more inclined than the iron 
shower of the same energy. Therefore, the EM compo- 
nent in the proton shower will attenuate more while reach- 
ing the ground from the shower maximum and it turns out 



A. YusHKOv et al. Precise determination of muon and EM shower contents 



18.50-18.60, 0-65° 




°'-e5' I 



u 





lg(E/eV) 



lg(E/eV) 



Figure 1: Top: ratio of ground plane signals at 1000 m 
in water Cherenkov detectors S'^/S'om vs vertical depth of 
shower maximum Xj^j^x ™^ 10^*'^ — lO^^'^ eV energy range 
for QGSJET II. Black line is the fit in the form (1). Bottom: 
means and RMS of distributions of relative difference be- 
tween MC simulated muon signals S^'^ and muon signals 
derived from the fit (1) 5fj*, calculated with the unique set 
of parameters for all energy bins: A = 538, h = —0.25, 
a — —0.22. Protons — red squares, iron — blue crosses. 



Figure 2: Top: EM ground signals at 1000 m in water 
Cherenkov detectors vs cos(6') in 10^^'^ — 10^^'^ eV energy 
range and 9 — 18° — 63° zenith angle range for QGSJET II. 
Black line is the fit in the form (3). Bottom: means and 
RMS of distributions of relative difference between MC 
simulated EM signals S^!^ and EM signals derived from 
the fit (3) S'^4, calculated with the unique set of parameters 
for all energy bins: 5°^^ = 2.53, co = —3, ci = 0.96, 
A = 0.012. Protons — red squares, iron — blue crosses. 
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that S'f^/S'P„ ratio becomes almost equal to the Sf/S^ total signal in water Cherenkov detectors: 
one, that allows to state a new shower universality prop- 
erty: the ratio of the muon signal to the EM signal S'p/5cm 
is the same for all showers, reaching the maximum at the 
same vertical depth X"^^^^, independently on the primary 
particle nature, primary energy and incident zenith angle 
(for the energy and angular ranges considered here). This 
property is illustrated in Fig. 1, where the dependence 
of Sfj,/Scm on X'^^^ for p and Fe primaries is shown in 
lg(£'/eV) = 18.5 — 18.6 energy bin. The functional de- 
pendence between X^^^ and S/^i/S cm turns out to be very 
simple and quasi-universal for all energies and primaries. 
The following function 



(2) 



-'^max — ^{Sfj./Sc 



(1) 



has been used to fit the data in 15 energy bins 
Alg(ii^/eV)=0.1 and the fit parameters have been found 
to be stable across the entire energy range. Using the 
functional dependence of S'^/S'om on X^-^^^ and Siqoq = 
5cm + S^ one easily gets the equation, which allows to ob- 
tain the muon signal from vertical depth of maximum and 



l + l/{{X^,JAy"'-a) 

We calculated the difference between the Monte-Carlo 
(MC) simulated muon signal Sjf''^ and the muon signal ob- 
tained from the fit S^^. In Fig. 1 we plot the behaviour of 
the mean and RMS values of these distributions for vari- 
ous energies, obtained with the unique set of fit parameters 
A — 538, b — —0.25 and a — —0.22, representing the av- 
erages over 15 A lg(_E/eV)=0. 1 energy bins. It is seen that 
the estimates of muon signals are unbiased with less than 
1% deviation of the mean reconstructed muon signal from 
the MC one for all primaries and the RMS values are small: 
8% for protons and around 5% for oxygen and iron (though 
we don't show results for oxygen, we use oxygen showers 
together with proton and iron ones to perform fits). 

In case of EPOS 1 .99 the same universality holds and the 
fit in the form (1) also provides good description of the 
simulated data, but, as expected, the coefficients of the fit 
are different from those for QGSJET II. 
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EM halo from muon decays and larger number of muons in 
iron showers brings to a larger EM halo signal. 

Looking at the showers at different zenith angles one sam- 
ples longitudinal showers profiles, for this reason it is nat- 
ural to try to describe the dependence of the EM signal on 
cos(6') with Gaisser-Hillas type function, using cos{6) as 
variable instead of X^ax: 
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cos(6') — Co 
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X exp 



ci - Co 
ci — cos(6') 
A 



(3) 



where a — {ci — co)/A; S^^^ (signal at maximum), co, 
ci (cosine of angle at which Scin=S^ni) and A are fit 
parameters. The fit parameters S^!^-^ and ci change by 
less than 10% and 3% correspondingly across the entire 
range of energies (when one makes fits in 15 energy bins 
Alg(^/eV)=0.1 from lg(£;/eV) = 18.5 to lg(£;/eV) = 
20.0), while cq changes quite chaotically from to —20 
(this causes A to change also). We have found that fixing 
Co (similarly to [9]) to any negative value within this range, 
we obtain a good universal fit and A changes in this case by 
less than 15%. Finally, we used the following average val- 
ues (except for cq that was fixed to —3) of the coefficients 
Sl!^-^ = 2.53, Co = -3, ci = 0.96, A = 0.012. The results 
of the fit and the difference between the MC simulated EM 
signal S^ and the EM signal obtained from the fit S^^ 
are shown in Fig. 2. The accuracy of the EM signal repro- 
duction for all energy bins is such that one gets an unbiased 
estimate of 5cm with RMS below 15% for proton and 13% 
for iron showers. 

Our calculations demonstrate that the universality of EM 
signal dependence on zenith angle holds true also in case 
of EPOS 1.99. 



Figure 3: Top: EM muon halo fraction S'cm,haio of total 
EM signal Scm vs zenith angle. Bottom: S/^JScm depen- 
dence on the zenith angle. lg{E/eV) = 18.50 - 19.50. 

2 Showers at the same zenith angles 

Another universality property follows from the study of 
showers arriving at the same zenith angles. In this case 
the average iron shower has to cross larger slant distance 
from X,„ax to the ground with respect to the average pro- 
ton shower and this almost equalizes EM signals for both 
primaries at the observation level in a wide range of zenith 
angles. For the signal at 1000 meters in the Cherenkov 
water detectors notable discrepancies between p and Fe 
EM showers components are observed for nearly vertical 
showers (9 < 18°, cos^(6') > 0.9) and very inclined ones 
(9 > 63°, 008^(61) < 0.2). In the first case the path from 
-'^max to the ground for p and Fe showers is almost the 
same. For inclined showers the difference is caused by the 



3 S^/Scm universaUty in respect to interac- 
tion models for ^ > 45° 



Phenomenologically the angular region 45° — 65° is of in- 
terest since with increase of the zenith angle the EM com- 
ponent produced mostly in vr*' decays at the initial EAS 
development stages is largely absorbed in the atmosphere 
and EM halo from muon decays starts to play a remarkable 
role (Fig. 3). One expects in this case that the behavior of 
the Sfj^/Scm ratio should become less sensitive to the prop- 
erties of the interaction models since with increase of the 
angle it more and more reflects the equilibrium state be- 
tween muons and EM halo from muons decays and inter- 
actions. To illustrate quantitatively this process let us write 
the S'p/S'om ratio for QGSJET II as 



CQGS / CQGS 
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here S'^n'^^aio i^ the EM halo signal from muons, SgJJfpure 
is EM signal from everything else except muons. Then for 
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EPOS 1.99 one gets 
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the scaling factors between muon and EM signals of the 
models and we have taken into account that 5*0111, halo oc 5^ 

In these notations one gets 



andsoM = 5-«i^aio/5. 

OEPOS / CEPOS 

^u / *^cm 



cm. halo' 



cQGS 



/5. 



QGS 

cm 



/qQGS /qQGS 

/^Wcm.halo/'-'' 



(4) 



cm, purcy 



One can see from Eq. (4) that with the increase of the zenith 
angle and hence of 5'^m,\aio/'5'?m,''purc the difference be- 
tween models in S^/Sem is decreasing as shown in the bot- 
tom panel of Fig. 3. Let us note that from the approximate 
equality of Sf^/Scm ratios for different models it follows 
that for any primary nucleus (p, O, Fe etc.) the following 
equality holds 



CQGS 

•^1000 



^EPOS 
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CEPOS 
cm 
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This, in turn, means that Eq. (2) in this angular range pro- 
vides an almost model-independent estimate of the muon 
signal. In fact, since the muon signal scales in the same 
way as the total signal, if one applies e.g. Eq. (2) with 
fit parameters for QGSJET II to the data simulated with 
EPOS 1.99, the total signal of EPOS 1.99 will give correct 
normalization for the muon signal. The difference between 
models in X^^^ and in the functional dependence on ^^^x 
will play only a minor role. As it will be demonstrated else- 
where [12] the muon signal for EPOS 1.99 can be retrieved 
with the use of the QGSJET II fit parameters with accuracy 
of 3-5%. 



advantage of 100% ground array duty cycle with respect to 
10% one of the fluorescence telescopes. On the other hand 
for angles below 45 degrees the difference in behaviour of 
S^/Scm between models should be large enough so that 
with simultaneous knowledge of S*^, S'em and Xmax one 
could be able to check the predictions of hadronic models 
quite easily using (1) and comparing e.g. the parameterized 
Xmax with the measured one. 

Finally, we would like to dwell on the problem of muon 
excess in the real data compared to predictions of the in- 
teraction models [9-11]. Since the muon content of EAS 
is highly model-dependent and the UHECR mass compo- 
sition is still unknown, this muon excess can be expressed 
only in terms of a relative excess with respect to the predic- 
tion of a given hadronic interaction model for a given pri- 
mary like (real signal)/(MC signal for protons). In [12] we 
show that for the zenith angles above 45 degrees it is pos- 
sible to get muon shower content from real data in almost 
interaction model independent way (see application to the 
data of the Pierre Auger Observatory in [13]). This, in turn, 
provides us with EM shower content which is weakly sen- 
sitive to the mass of the primary particle and hence allows 
to find the absolute scaling factor (real signal protons)/(MC 
signal protons). Performing such scaling on the one hand 
will diminish the difference in predicted muon signals be- 
tween different interaction models, and on the other hand 
will open the possibility to use the muon signal for mass 
composition studies. 
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Conclusions 

We have presented two new EAS universality properties 
providing two independent ways to access EM and muon 
shower contents. We have shown that these properties can 
be described with simple parametrizations which are valid 
in wide energy and zenith angle ranges, and are indepen- 
dent on the primary particle nature. We believe that these 
universality properties can be used in hybrid experiments 
for mass composition studies, for primary and missing en- 
ergy estimates and for tests of hadronic interaction models. 
One of the possible strategies lies in the simultaneous ap- 
plication of both universality properties to the data. It is 
clear that parametrizations (2) and (3) will give consistent 
estimates of muon and EM shower contents only in case 
of correct description of the hadronic interaction properties 
by the particular model. Another interesting strategy can be 
pursued in hybrid experiments equipped with muon detec- 
tors. For zenith angles above 45 degrees where the EM halo 
plays an important role, this universality property can be 
used for the determination of the depth of the shower max- 
imum in almost interaction model independent way taking 
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